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We construct and analyze a microscopic model for insulating rock salt ordered double perovskites, with the 
chemical formula A2BB'C>6, where the magnetic ion B' has a 4d 2 or 5d 2 electronic configuration and forms 
a face centered cubic (fee) lattice. For these B' ions, the combination of the triply-degenerate antisymmetric 
two-electron orbital states and strong spin-orbit coupling forms local quintuplets with an effective spin moment 
j — 2. Moreover, due to strongly orbital-dependent exchange, the effective spins have substantial biquadratic 
and bicubic interactions (fourth and sixth order in the spins, respectively). This leads, at the mean field level, to 
a rich ground state phase diagram which includes seven different phases: a uniform ferromagnetic phase with 
an ordering wavevector p = and uniform magnetization along [111] direction, four two-sublattice phases 
with an ordering wavevector p = 27r(001) and two four-sublattice antiferromagnetic phases. Amongst the two- 
sublattice phases there is a quadrupolar ordered phase which preserves time reversal symmetry. Extending the 
mean field theory to finite temperatures, we find ten different magnetization processes with different magnetic 
thermal transitions. In particular, we find that thermal fluctuations stabilize the two-sublattice quadrupolar 
ordered phase in a large portion of phase diagram. Existing and possible future experiments are discussed in 
light of these theoretical predictions. 

PACS numbers: 71.70.Ej,71.70.Gm,75.10.-b 



I. INTRODUCTION 

The combination of strong electron correlation and strong 
spin-orbit coupling (SOC) is relatively unexplored theoreti- 
cally. It arises naturally in a broad family of magnetic Mott 
insulating systems in which the three fold degenerate t^g or- 
bitals are partially filled. In these systems, the t2 g orbital de- 
generacy is protected by cubic lattice symmetry and the crys- 
tal field splitting is large enough so that e g orbitals are not 
occupied. Unlike for e g states, SOC is unquenched for ti g or- 
bitals and splits the one electron levels into an upper j = 1/2 
doublet and a lower j — 3/2 quadruplet™ 

In this category, many Ir-based magnets have been studied 
both theoretically and experimentally!^ Here, the magnetic 
ion Ir 4+ has a d 5 electron configuration with five electrons 
residing on the t2 g orbitals and the effective j = 1/2 descrip- 
tion can be adopted. Some very exotic states like a quantum 
spin liquid in Na^rsO^^and a topological Mott insulator in 
A2lr20'p'^'have been observed and proposed in the strong 
coupling and intermediate coupling regime, respectively. 

Moving beyond iridates, in our recent workj^l we have 
studied the magnetic properties of a series of compounds 
called ordered double perovskites with the chemical formula 
A2BB'06.^^' We considered double perovskites in which 
the magnetic ions B' (e.g. Re 6+ , Os 7+ , Mo 5+ ) have a d 1 
electron configuration with one electron residing on the ti g 
orbitals and hence form j = 3/2 local moments. In this anal- 
ysis, we found several exotic phases including a novel ferro- 
magnetic state driven primarily by orbital interaction, an anti- 
ferromagnetic state with strong octupolar order and a spin ne- 
matic state, and furthermore a quantum spin liquid state pos- 



tulated in a region of the phase diagram. 

To round out the list of of magnetic systems with par- 
tially filled t2 g orbitals, we must consider the d 2 , d 3 and d A 
cases. For a d 3 electron configuration, the three electrons fill 
all the three single-electron t 2g orbitals, forming an antisym- 
metric orbital wavefunction. The orbital degree of freedom is 
completely quenched. The system is described by spin-only 
Hamiltonian with spin S = 3/2. Since it has a large spin, 
one may expect it to behave rather classically.-^For a d 4 elec- 
tron configuration, when the SOC dominates over the Hund's 
coupling the four electrons completely fill the lower j = 3/2 
quadruplets and there is no local moment description at lowest 
order of approximation. When the Hund's coupling dominates 
over SOC, the four t2 g electrons form a total spin S = 1 but 
still have a three fold orbital degeneracy. The effective SOC 
further lifts the spin-orbital degeneracy completely and also 
favors a trivial j ' = local state. So the only non-trivial case 
left is the d 2 configuration, with two electrons filling the t2 g 
orbitals. 

In this paper, we consider this valence state in the context 
of double perovskites mentioned above, specifically extend- 
ing the theory of Ref. 9 to the case of a B' ion with a 4d 2 
or 5d 2 electron configurationPED Unlike for the 4c? 1 or 5c? 1 
system with one electron per site, the local Coulomb inter- 
action plays an important role in determining the the local 
spin and orbital structures for the Ad 2 or 5d 2 systems. For 
the spin sector, the first Hund's rule requires a symmetrized 
spin wavefunction, favoring a total spin 5 = 1. For the orbital 
sector, the orbital electron wavefunction should be antisym- 
metrized, composed of two single-electron t2 g orbitals. These 
three antisymmetrized two-electron states act as an effective 
I = 1 total orbital angular moment. The strong SOC com- 
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bines the total spin 5=1 with the total orbital moment I = 1 
and leads to an effective total angular momentum j = 2 de- 
scription of the system. Similar to the j = 3/2 case studied 
in Ref. 9 the orbitally-dependent nature of the spin interac- 
tions lead to an interesting microscopic Hamiltonian which 
contains significant biquadratic (fourth order in spin opera- 
tors) and triquadratic (sixth order in spin operators) spin-spin 
interactions. These unusual interactions may be understood 
as couplings between the local magnetic quadrupole and oc- 
tupole moments. Therefore, an analysis of the microscopic 
Hamiltonian naturally leads to a rich structure of magnetic 
multipolar orders. 

The results of a mean-field analysis are summarized in 
Fig. [T] and Table [I] There are seven total ground state phases 
which appear, including notably a broad region of time- 
reversal invariant but quadrupolar ordered (spin nematic) 
ground state with a two-sublattice structure of the orbital con- 
figuration. This is described by the quadrupole tensor opera- 
tors, 



qf = mf - a?) 2 - (i?) 2 ]/Vs, 

Qf- y2 = (jf) 2 - (j!) 2 - 



(i) 

(2) 



In the quadrupolar phase, (Q^ ~ v ) = ±q' alternates sign on 
the two sublattices (see Eqs. {73])). In this phase, the time 
reversal symmetry is unbroken and there is no magnetic dipo- 
lar and octupolar order. Such a local spin nematic ground 
state is particular to integer spin systems and prohibited for 
half-integer spins. A similar spin nematic ground state has 
also bee n prop osed theoretically for a spin 5 = 1 material 
NiGa2S4p^"24lbut an experimental confirmation of this phase 
in NiGa2S4 is still lacking. Apart from the quadrupolar or- 
dered phase, the other ground states are magnetic, and com- 
prise both ferromagnetic and antiferromagnetic states with en- 
larged unit cells up to quadruple the size of the ideal one. The 
extension of the mean field analysis to T > is described 
in Fig. [2] Amongst various finite temperature phases, the 
quadrupolar ordered one is prominent, providing more oppor- 
tunity for its experimental discovery in real material. 

The remainder of the paper is organized as follows. In 
Sec. |ll| we first explain the on-site spin and orbital physics 
which leads to an effective j — 2 local moments on each B' 
magnetic ion. As we did in previous workp we introduce a 
microscopic Hamiltonian which includes three interactions: 
nearest neighbor (NN) antiferromagnetic (AFM) exchange, J, 
NN ferromagnetic (FM) exchange, J', and electric quadrupo- 
lar interaction, V. Due to strong SOC, we project these in- 
teractions, which are written in the separate spin and orbital 
spaces, down to the j — 2 manifold. This leads to a spatially 
anisotropic and very non-Heisenberg-like Hamiltonian, which 
contains many terms beyond the usual quadratic exchange. In 
Sec. [Hi] we find the mean field ground state phase diagram 
of the model Hamiltonian and analyze the properties of each 
ground state. Some insight into the general results is given 
by considering some simple limits, including both strong easy 
plane and easy axis anisotropy, and an orbital interaction only 
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FIG. 1. (Color online) Ground state phase diagram of the model 
Hamiltonian in Eq. {24). Here V is the electric quadrupole interac- 
tion, and J' is the nearest-neighbor ferromagnetic exchange between 
orthogonal orbitals; both are normalized to a unit value, J = 1, for 
the nearest-neighbor antiferromagnetic coupling. The names for the 
phases are defined in Table]!] and in Sec. |III| 
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TABLE I. Phases in mean field theory. The first column gives the 
name of the phase, as shown in Fig. [T] In the second column, m 
denotes the direction (up to symmetries) of the uniform magnetiza- 
tion; a "0" is shown if this is vanishing. Similarly, in the third col- 
umn, n denotes the direction of the staggered magnetization. For 
the quadrupolar phase, in which "0" appears in both the second and 
third columns, time-reversal symmetry is unbroken. The fourth col- 
umn specifies the number of sites in the magnetic unit cell. The fifth 
column gives the number of site in the unit cell for time-reversal in- 
variant observables, i.e. the quadrupolar tensor. 



for the complete model Hamiltonian and find the zero temper- 
ature ground state phase diagram which is depicted in Fig. [T] 
In Sec. 



IV 



model (Sec. IIIBi. We find that, besides a uniform orbitally 



ordered Finally, in Sec. Ill C we carry out a mean field study 



we extend the mean field theory to T > 0, and 
identify the structure of the magnetic order for each phase. We 
find ten different finite temperature magnetization processes 
which correspond to the ten overlapping regions between the 
ground state phases and the shaded areas in Fig. [2] We also 
discuss the different finite temperature phase transitions asso- 
ciated with magnetic dipolar and quadrupolar orders. Tran- 
sition temperatures are extracted from Landau theory. In a 
broad region of the parameter space, we find the quadrupo- 
lar phase occurs at the intermediate temperatures between the 
high temperature paramagnetic phase and various low temper- 
ature magnetic ordered phases. 

Finally in Sec. [V] we compare our theoretical predictions 
with current experimental findings and suggest further direc- 
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tions for theory and experiment. 



H. MODEL 

A. Spin-orbit interaction and electron orbitals 

The magnetic ions B' (Re 5+ , Os 6+ ) in the relevant ordered 
double perovskites (Ba 2 CaOsO e , La 2 LiRe0 6 , Ba 2 YRe0 6 ;P 
all have a Ad 2 or 5cP electron configuration with two electrons 
on the triply degenerate t 2g multiplets. Because of the electron 
interaction, the local spin orbital state is quite different from 
the case of d 1 electron configuration where the single electron 
state is enough to describe the local physics. Considering the 
dominance of the crystal field splitting over the SOC, we now 
fill the three t 2g orbitals with these two electrons before in- 
cluding the effect of SOC. To respect the first Hund's rule, the 
total spin for the two electrons is S = 1. For the orbital sector, 
there are three degenerate antisymmetric two-electron states, 
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(\xy) 1 \xz) 2 - \xy) 2 \xz) 1 ) (3) 
{\xy) 1 \yz) 2 -\xy) 2 \yz) 1 ) (4) 
(\xz)i\yz) 2 - \xz) 2 \yz)i) , (5) 



in which, the subindex ("1" and "2") labels the electron. 
Therefore, there are totally nine-fold spin-orbital degenera- 
cies. The presence of SOC will lift some of the degenera- 
cies. Following the spirit of degenerate perturbation theory, 
we project the SOC onto the triplet subspace spanned Eq. d3), 



%so = — A 1 • S 



(6) 



in which, the total angular momentum quantum number of 
these operators are I = 1,5*= 1. The effective orbital angular 
momentum 1 comes from the projection of the total orbital 
angular momentum L = Li + L 2 onto the tripets in Eq. ((5), 



V L V = -1 



(7) 



Here V Q = X)a=x yz I ^) (M ^ s tne projection operator to the 
triplet orbital subspace. 

The reduced SOC in Eq. (|6]l favors a local j = 2 (j — 
1 + S ) over other higher energy states j = 0, 1 by an energy 
separation 0(A). In the materials we are considering, the SOC 
A is a very large energy scale (some fraction of an eV). 

In general, cubic symmetry allows the presence of an on- 
site cubic anisotropy term, (j x ) 4 + {j y ) i + (j 2 ) 4 , which lifts 
the degeneracy of the five j = 2 states. However, we ex- 
pect this splitting to be rather small, and provided it is smaller 
than the typical exchange coupling between spins, the j = 2 
description should be a good approximation. Microscopi- 
cally the cubic anisotropy comes from the 4th order effect of 
the SOC and pair hopping (between different orbitals on the 
same ion) terms J p which excite the electrons into the e g or- 
bitals. The magnitude of the cubic anisotropy should be of 
- C(A 4 /A 3 , J 4 /A 3 ) (with A the crystal field splitting be- 
tween e g and t 2g levels). This is certainly a much smaller 



energy scale compared to SOC, and likely small compared to 
exchange. In any case, we will neglect it in the following. 

In the strong SOC limit, every local operator should be pro- 
jected onto the local subspace spanned by five j = 2 states. In 
particular, 



v 2 sv 2 = l -i 
v 2 \v 2 = l - j 



(8) 
(9) 



Here V 2 is the projection operator into the local j = 2 states. 
In addition, one can find the local magnetic moment is given 
by 



M = V 2 (2S - 1) V 2 = *j 



(10) 



hence, the magnitude of the local magnetic moment is found 
to be y/6/2/j, B sy 1.25/zb- 



B. Exchange interactions and electric quadrupolar interaction 

In this subsection, we introduce the interactions between 
the local moments. From previous workj^we will need to con- 
sider the nearest neighbor (NN) antiferromagnetic (AFM) ex- 
change, NN FM exchange and NN electric quadrupolar inter- 
actions, and these interactions are highly anisotropic in both 
the position space and spin space. For example, in the XY 
plane, only electrons on xy orbital can virtually transfer from 
one site to another via the intermediate oxygen p orbitals. 
Thus, one finds that the NN AFM exchange is written as 



n/XY 



J 



E [a 

(tj)exy 



1 

— n. 
4 



t.xy' ' ;..; // | 



(ID 



where the sum is over nearest neighbor sites in the XY planes, 
and the correponding terms for YZ and XZ planes can be ob- 
tained by the obvious cubic permutation. One should note that 
the operators Sj ^ and rii <xy denote the electron spin residing 
on the single-electron xy orbital and orbital occupation num- 
ber for the single-electron xy orbital at site i, respectively. 
To connect these single electron operators to the two-electron 
operator which acts on the two-electron orbitals in Eq. d5), we 
have the following relations 

n iiXV = n;,x + n itY = {It) 2 (12) 

S w = |Kx+n, Y )=y(i 1 f ■ (B) 

Here n^x (or n^y) denotes the occupation number for |X) (or 
|Y)) of the two-electron orbital states at site i, and Si is the 
total spin 5=1 for the two electrons. The physical meaning 
of Eq. ( fT3~| is apparent. The electron occupation number on 
the single-electron orbital xy can be nonvanishing only when 
the two-electron orbital state |X) or | Y) is occupied by the two 
electrons. 

Throughout this paper, we use the subindices (i, xy) to de- 
note the site and single electron orbitals, subindex X to denote 
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the two-electron orbitals, superindex (/i = x,y,z) to denote 
the spin component, and capital letters (XY, XZ, YZ) to de- 
note the planes. With these definitions, we note the double 
occupancy condition at each site, which defines the Mott in- 
sulating phase, becomes, 



+ rii 



(14) 



in terms of the two-electron operators, the above equation is 
equivalent to 

rii,x + ni,Y + n hZ = 1 • (15) 
Moreover, from Eq. (13) , orbitally-resolved spins satisfy 



jection for Si >x -y and ni tXy , 



ox _ 1 -x | r -x\3 
i,xy 10 19 * ' 



\2 Jl 12 

i.xy 12 12 ' 



— i' x 4- — (i^ 3 

i,xy 12 12 ' 



(20) 
(21) 
(22) 
(23) 



in which, O = V-zOV?.. After the projection, the minimal 
Hamiltonian that we will study in this paper is 



ft — %A¥M + %FM + %quad- 



(24) 



^i,xy ~t~ ^i,xz &i,yz — • (16) 



The second interaction to include is the NN FM exchange 
interaction. FM exchange comes about when the orthogonal 
p orbitals on a single oxygen ion are involved in the exchange 
path. Directly using the results from Ref. 9 and the relation 
in Eq. ( fT3] >, one can immediately write down this interaction. 
Again, for two sites i,j in the XY plane, this FM exhange is 
given as 



njXY 
n FM,ij 



J \^i,xy ' (&j,yz + Sj )XZ ) I (l f> j) 

3J' 



2 

,r_ 

4 

3J' 



+ j) 



The third interaction to include is the electric quadrupolar 
interaction. This is obtained by evaluating the Coulomb inter- 
action in different orbital occupations. As the AFM and FM 
exchange, we also take results from previous work to write 
down this interaction. In XY plane, we obtain the electric 
quadrupolar interaction as 



H, 



XY 

quad, ij 



W 

X 

9V 

~~^~ n i,xy n j : xy 
W 

9V 



W'i^yz) {j^j.xz Tlj.yz ) 



+ "-jr(lt) 2 W 2 



(18) 



The minimal Hamiltonian for the cubic system contains all 
three of these interactions in addition to the onsite SOC, 



H = H 



AFM 



FM 



^quad + ft 

SO 



(19) 



Since we are interested in the limit of strong SOI, we have 
to project the minimal Hamiltonian H onto the five j = 2 
states at every site. As an illustration, we write down the pro- 



Ill. MEAN-FIELD GROUND STATES 

In this section, we study the zero temperature phase dia- 
gram of the model Hamiltonian in Eq. d24]i. Ultimately, in 



Sec. IIIC we will do this by mean field theory, or equiva- 
lently, variationally searching for direct product states which 
minimize the expectation of W. Before reporting these results, 
however, we discuss some simple limits in which the behav- 
ior can be understood more intuitively. First, in Sec. Ill A 



we impose a strong single-ion uniaxial anisotropy, which re- 
moves the orbital degeneracy renders the problem trivially sol- 
uble with singlet, ferromagnetic and antiferromagnetic ground 
states. Second, we consider a pure orbital model, in which 
only the electric quadrupole interaction V is included. This 
gives the two sublattice quadrupolar state described in the in- 



traduction. Finally, in Sec. Ill C we report the results of a full 



mean-field calculation including all couplings and making no 
further approximations. 



A. Uniaxial anisotropy 

As we did in Ref. |£l we first consider the ground state of 
this Hamiltonian in the presence of strong easy-plane or easy- 
axis anisotropies. The strong easy-plane anisotropy (on XY 
plane) is a trivial limit and is modeled by ^ D(jf) 2 with a 
positive D. When D is quite large (compared to exchange 
coupling and electric quadrupolar interaction), the spin state 
on every site is pinned to \ j z = 0), which is a rather triv- 
ial uniform state with an ordering wavevector p = 0. The 
strong easy-axis anisotropy (along z direction) is less trivial 
and is modeled by the same Hamiltonian but with a negative 
D. Large \D\ favors either \j z — 2) or \j z = —2) to be 
occupied. After projecting the Hamiltonian in Eq. ( |2"4"| ) onto 
this two states, the electric quadrupolar interaction is com- 
pletely quenched and the resulting effective Hamiltonian is 
a trivial Ising Hamiltonian. One can readily find that, when 
J' > 5 J/38 the ground state is a ferromagnetic state with 
an ordering wavevector p = 0, and when J' < 5 J/38 the 
ground state is an antiferromagnetic state with an ordering 
wavevector p = 27r(100) or 27r(010). One may postulate 
from these anisotropic case that, the ground state for the actual 
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cubic Hamiltonian may either have a uniform state (p = 0) or 
a two-sublattice state (p = 27r(001) and equivalent wavevec- 
tors). As we will see in the following sections, this guess is 
correct for a large portion of the parameter space, but we also 
find some interesting exceptions. 



B. The orbital Hamiltonian 

To understand the nature of the non-magnetic quadrupolar 
ground state of the model, it is sufficient to consider only the 
electric quadrupole interaction V, given in Eq. ( (31) . Note that 



Tli,xz and Tli^xy 



it involves only the three operators n$ . 
each site. Since these are all time-reversal invariant, it is ap- 
parent that they do not span the full space of j = 2 operators. 
Thus there must be additional constants of the motion, and 
the Hamiltonian can be separated into sectors corresponding 
to different irreducible representations (irreps) of the algebra 
of these operators. Indeed, one can show that the single-site 
Hilbert space decomposes into one two-dimensional irrep and 
three one-dimensional irreps. The two dimensional irrep is 
spanned by the two states 

1 :ii«=2)+u* = -2)) 



1 ' V2 
\d) = \f=0). 
In this subspace, the orbital operators become 



rii- 



1 

2y/t 
1 

2^3 ( 



where (?j are the Pauli matrices acting in the \u) , | d) 
the first of the three one dimensional irreps we have 



(25) 
(26) 

(27) 

(28) 

(29) 
space. In 









^i,xz 1 




(30) 









The other two one dimensional irreps may be obtained by per- 
muting these values. 

The ground state must consist of a single irrep on each site. 
Which irrep occurs should be determined by minimization of 
the energy. Consider the case in which each site has a one- 
dimensional irrep. Then we can specify the state of a site by a 
"Potts"-like variable s, = 1,2,3 specifying which of the three 
orbital numbers equals 1, i.e. 

with a = 1,2,3 corresponding to a = yz,xz,xy, respec- 
tively. Then the Hamiltonian for a bond in the XY plane be- 
comes 



XY 



V 



1 d irreps 



3 
9V 
" 16 



(1 



with YZ and XZ plane bond interactions obtained by permuta- 
tions. In this case the Hamiltonian is purely classical and thus 
the ground state can be exactly found by minimization. We 
find that a ferromagnetic ground state is preferred, with con- 
stant Si (there are thus three such degenerate ground states), 
which has an average energy of 65^/72 sw 0.903V^ per bond. 

The other natural choice to consider is when the two- 
dimensional irrep is chosen on each site. In this case the 
Hamiltonian can be written as 



XY 



2d irrep 



V 



25 „xx 

36^ 



1 



(33) 



where the dot product involves only the x and z components 
of the Pauli matrices, and we have dropped a term linear in the 
Pauli matrices which cancels when all three types of bonds are 
added. In this case the interactions for bonds in the other two 
planes are obtained by the rotations erf — > —\<j\ 



J- 2 u. 



with the upper (lower) sign chosen for the XZ (YZ) plane. 
The dot product is unchanged by this rotation. 

The Hamiltonian in this subspace is a type of Kugel- 
Khomskii model, similar to that studied in models of e g or- 
bitals. It is fully quantum, and thus cannot be solved ex- 
actly. However, within mean field theory we find the varia- 
tional ground state is simply the anti-ferromagnetic product 
state with erf = (— l) z , i.e. with alternating sign on adjacent 
XY planes. In this state the expectation of the Hamiltonian 
can be taken, and the energy is found to be E v = 173V^/216 rj 
0.80 IF per bond. Note that this is lower than the energy found 
for the one-dimensional irreps. This is an upper (variational) 
bound on the ground state energy in this sector, so we indeed 
expect the ground state to be in the two-dimensional irrep. 

Physically, the mean-field state describes an orbitally or- 
dered phase with a two sublattice structure with single-ion 
wavefunctions 

I A) = \\f = 2) + ±=\f = 0) + \\f = -2), (34) 

IB) = \\f = 2) - j=\j* = 0) + \\f = -2), (35) 

where A, B label the two different sublattices (planes with 
even/odd z). Remarkably, these states are invariant under time 
reversal. Since time reversal symmetry is unbroken, there is 
no magnetic order, 



0.) = o. 



(36) 



Therefore, this phase is a magnetic quadrupolar phase (or spin 
nematic phase) and the orbital configuration is given by 



« _ ,2 , n/3 2 ^3 2. 
(n A ) -(- + --, - ---,-) 



(fi B ) = G 



V3 2 
~6~' 3 



VS 2 
IT' 3 



(37) 
(38) 



3) , (32) 

Here, we defined n = (n yz , h xz , n xy ) for convenience. 
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C. The full cubic Hamiltonian 

Although both the anisotropic limit and pure orbital interac- 
tion support a two-sublattice ground state, it is still question- 
able that the cubic Hamiltonian will also behave likewise. In 
this section we report the results of a systematic investigation 
of the mean field ground states of the full Hamiltonian, allow- 
ing for large unit cells (we considered cells of up to 4 sites). 
We made no further assumptions and variationally minimized 
the energy with respect to an arbitrary wavefunction on every 
site of the unit cell. Finally we verified that each mean field 
ground state is stable within linear flavor wave theory!^ The 
mean field phase diagram is depicted in Fig. [T] and the key 
features of each phase is listed in Table [I] Within the flavor 
wave theory, all seven phases exhibit an energy gap. Next we 
describe each of the seven phases. 



1. Antiferromagnetic (AFM100) state 

For small J'/ J and V/J, the ground state is a typ- 
ical antiferromagnetic phase with an ordering wavevector 
p = 27r(001). States with the equivalent momenta p — 
27r(100), 27r(010) are of course degenerate. To be specific, 
we will take p = 27r(001) for all the two-sublattice phases in 
the following. The variational mean field ground state wave- 
function on the A and B sublattices has the following form: 



\A) = —(\f=2) + \f 
+ \{\f = l) + \f-- 

|B) = -^(|j*=2) + tf a 
-\{\f = l) + \f-- 



-2)) 

-1» 
-2)) 



1 



x 2 \f = 0) 



(39) 



x 2 \f = 0) 



(40) 



in which x is a real parameter, which is found by minimiz- 
ing the variatonal energy. Since under time-reversal, \j z — 
to) — > (— l) m \j z = m), these two states transform into one 
other under time reversal. Therefore, the magnetic dipolar and 
octupolar orders are anti-parallel on two sublattices. Because 
of the intrinsic strong SOC, the local moments are aligned by 
crystalline anisotropy, and in this state orient along the [100] 
(or equivalently, [010]) axis. A more precise description of 
the symmetry breaking of the phase is given by introducing 
the magnetic dipole and quadrupole moment operators. The 
magnetic dipole moment is antiferromagnetically ordered, 



(j A/B )=±m(l,0,0), 



(41) 



with m = \/2x + -\/3(l — 2x 2 ). The magnetic quadrupole 
tensor, however, is uniformly ordered, 



(Qf ) = 2\/3(3(n 

(Qf- y2 ) = m 



/..((/ 
ft 



2)=q 
: ... z j Q 



(42) 
(43) 



with q = 2V3{2x 2 



3/4) and q' = 2^/3(1 - 2x 2 ) + 3/2. 



The uniform quadrupolar order (or orbital configuration) 
can be understood to arise from the large NN AFM exchange 
J which favors time reversal pairs on two sublattices. As the 
ferromagnetic exchange and electric quadrupolar interaction 
increase, the orbital-orbital interaction will become important 
and the uniform orbital structure will break down. 



2. Uniform ferromagnetic (FM111) state 

With large J' / J and small V/J, the ferromagnetic ex- 
change dominates and favors a uniform ground state with the 
spin polarization aligned with [1 1 1] or other equivalent lattice 
directions. The mean-field ground state of this phase is a fully 
polarized spin eigenstate with quantization axis along [111], 
so that 



n/3 



(1,1,1), 



(44) 



with to = 2. And the three orbitals are equally populated, 



n, 



2 2 2 
l 3' 3' 3 j ' 



(45) 



(46) 



and the magnetic quadrupolar orders vanish, 

(Qf) = (Qf- y2 )=0- 

3. Two-sublattice ferromagnetic (FM110) state 



With large J' / J and V/J, we obtain a FM1 10 state which 
is the same phase proposed in Ref. 9 This state can be con- 
sidered as a compromise between the tendencies of J and J' 
to order the moments along the [100] and [111] axes. The 
competition between these two effects allows the orbital inter- 
action to stabilize a coexisting quadrupolar order. The result 
is that the orbital configuration has a two-sublattice structure, 
and the magnetic moments on the two sublattices are neither 
antiparallel or parallel. The ground state wavefunction of this 
phase is parametrized by two complex numbers, x\ and x%, 

1 ' S2|i* = -2)) + a;o|i* = 0> 



|A> = 



V2 
1 

71 
i 



{x 2 \f = 2) 
(xi\f 



= 1 



b* = -i» 



(47) 



B) = -^{-nf = 2) - x 2 \f = -2))+ x \f = 0) 



V2 

+ ^(-m\f = i) 



IX\ 



\f- = -1)) 



(48) 



where, xq = y 1 — \ x i\ 2 ~ |^2p and xY/z is complex conju- 
gation of Xx 2- From the wavefunction, we find that magnetic 
dipole and quadrupole moments have the following form: 

:(l,l,0)±m'^(l,-l,0) 



(Q 



Ua/b) 



A/B 



V2 

= q 

= ±q>. 



V2 



(49) 

(50) 
(51) 
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Here, the actual expression of m, m', q, q' in terms of X\ and 
x 2 is quite involved and not important for the purpose of pre- 
sentation. So we will not write them out explicitly. Similar 
omissions are made in later sections. 

From Eq. ( |5"Tj ), the total magnetic moment is along [110] di- 
rection, and the staggered magnetic dipole moment is perpen- 
dicular to the total magnetic dipole moment and along [110] 
direction. The orbital configurations of the two sublattices are 
similar to that in the quadrupolar phase (see Eq. ([38))). 



4. Intermediate ferromagnetic ("* ") state 

Between the FM111 and FM110 phases, we find an in- 
termediate state, which we denote This state also has 
a two-sublattice structure with the ordering wavevector p = 
27r(001). We find the ground state wavefunction of this phase 
is parametrized by four complex numbers, Xx ) X2,X-i and 

X-2, 

| A) = x 2 \f = 2) + x x \f = 1) + x \f = 0) 

+ x^\f = -1) + x- 2 \f = -2) (52) 

|B) = -Sab* = 2) - ixi\f = 1) + x \f = 0) 

+ ixZ^\f = -1) - = -2) (53) 



with x = y/l - \xi\ 2 - \x 2 \ 2 - \x-i\ 2 - \x- 2 \ 2 . The mag- 
netic order in this phase interpolates between that of the 
FM111 and FM110 phases. The uniform magnetization ori- 
ents along an axis between the [111] and [110] directions. 
Like in FM1 10 phase, the staggered magnetic dipole moment 
in phase "*" is along [110] direction. The orbital configura- 
tion has the same two-sublattice structure as in FM1 10 phase. 
The corresponding order parameters are 



Oa/b) = (mi, mi, m 2 ) ± m' -^=(1,-1,0), (54) 



(QT/b) = 9, 

(Qir/)=±<i'- 



(55) 
(56) 



parametrized by two complex number, x% and x 2 , 

I A) = ^2\f = 2)+x- 2 \f = -2))+ x \f = 0) 
+ -^( Xl \f = l)+xl\f = -!)) (57) 

IB) = = 2)+x- 2 \f = -2)) + x \j* = 0) 

+ -^(-x 1 \f = l)-x- 1 \f = -l)) (58) 

|C) = ^(-^\f = 2 > - = -2)) + x \f = 0) 



+ -^(-ix x \f = 1) + ixj\f = -1)) 



(59) 



ID) = ^(-*2\f = 2 > - = -2)) + x Q \f = 0) 



+ —( tXl \f = 1) - txl \f = -1)) 



(60) 



with x = y^l — \xi\ 2 — \x 2 \ 2 . The order parameters are: 

Oa/b) = ±m(«i,u 2 ,0), (61) 
Oc/d) = ±m(-u 2 , ui, 0), (62) 

(Qa/b/c/d) = ^ (63) 

2 2 „2 2 

(64) 



(Ql/B y > = (Q x c/ u v > = ±?'- 

It is easy to see from the above order parameters that the 
orbital configuration still has a two-sublattice structure as de- 
manded by the orbital-orbital interaction while the magnetic 
dipolar order has a four-sublattice antiferromagnetic structure. 
One can think of this state as breaking time reversal symmetry 
on top of a two-sublattice orbitally ordered state by develop- 
ing antiferromagnetic orders within each sublattice. States on 
sublattice A and B form a time-reversal pair and states on sub- 
lattice C and D form another time-reversal pair. But it is not 
a conventional antiferromagnetic state. In fact, in this phase, 
the magnetic dipolar moment of A or B sublattice is nearly 
perpendicular to the magnetic dipolar moment of C or D sub- 
lattice. Consistent with strong magnetic anisotropy, the stag- 
gered magnetizations (jA-jB+jc-jD) and (jA-jB-jc+jD) 
are oriented in a direction very close to [110] lattice direction. 



5. Four-sublattice antiferromagnetic ( "A ") state 

In the region of a small J' / J and an intermediate V/J, 
the FM exchange interaction between time-reversally odd mo- 
ments has negligible effec, while the AFM exchange and elec- 
tric quadrupolar coupling are somewhat balanced. The keen 
competition between these two interactions induces an in- 
teresting intermediate phase: a four-sublattice antiferromag- 
netic phase which we denote" A". The magnetic unit cell 
is the elementary tetrahedron of the fee lattice. Similar to 
FM110 phase, the ground state wavefunction is found to be 



6. Four-sublattice antiferromagnetic ("A") state 

In the intermediate J' / J and the small V/J regime, we 
find another four-sublattice antiferromagnetic phase, the "A" 
state. In this regime, the electric quadrupolar interaction may 
be neglected, and we can understand the state as arising due 
to the competition between FM and AFM exchange interac- 
tions. In the AFM100 phase, every site has eight NN AFM 
neighbors and four NN FM neighbors. In the FM111 phase, 
every site has zero NN AFM neighbors and twelve NN FM 
neighbors. In the "A", we find that, for every site there is four 
NN AFM neighbors which is an intermediate case compared 
to FM1 1 1 and AFM100 phase. The ground state wavefunction 
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of this "A" phase is parametrized by two complex numbers, 
x\ and X2, 

|A) = fyvti* = 2) + xi\f = -2)) + x \f = 0) 

+ ^(x 1 \f = l)+xi\f = -1)) (65) 
IB) = ^2\f = 2)+x^\f = -2))+x \f = 0) 

+ ^{-xi\j' = l)-Wi\j' = -l)) (66) 
|C) = = 2) + x 2 \f = -2)) + x \f = 0) 

+ i(3T|j* = l>+s 1 |j* = -l>) (67) 
ID) = ^ W = 2) + x 2 \f = -2)) + .x |.f = 0) 

+ ±(-xj\f = 1) - Xl \f = -1)) (68) 



8. T = Transitions 



At the mean field level, the nature of the transitions between 
different phases can be understood from the wavefunctions. If 
the wavefunction of one phase can be continuously tuned to 
that of the neighboring phase, then the phase transition be- 
tween these two phases may be continuous. Otherwise, it 
is first order. We find the possible continuous transitions are 
FM111-"*", "*"-FM110, FMllO-Quadrupolar, Quadrupolar- 
"A" and "A"-AFM100. The remaining transitions are first 
order at mean field level. We do not discuss fluctuation ef- 
fects here, which would be required for a full understanding 
of the transitions. 



IV. T > PHASES 



with so = \J\ — |a;i | 2 — I22I 2 . As one can see from the 
wavefunction, states on A and B sublattices form a time rever- 
sal pair and states on C and D sublattices form another time 
reversal pair. The magnetic dipolar and quadrupolar orders of 
each sublattice are found to have the following relation, 



Oa/b) = ±m(ui,tt2,0), 
Oc/d) = ±m(ux, -i«2,0), 

Wa/b/c/d) = 1i 

(Q " 



x -y 
A/B/C/D- 



= q 



(69) 
(70) 

(71) 
(72) 



Although both "A" and "A" states are four-sublattice 
states, they are actually different phases. For instance, un- 
like in the phase "A" discussed in last section, the quadrupole 
moments in the "A" phase are uniform. Intuitively, this is 
because the orbital-orbital interaction is not large enough to 
induce a two-sublattice structure in this regime. 



7. Two-sublattice quadrupole (spin nematic) state 

With an intermediate J'/ J and a large V/J, the orbital- 
orbital interactions dominate over all other interactions in the 
Hamiltonian. Therefore the ground state reduces to quadrupo- 
lar phase found in Sec. |IIIB| The full set of order parameters 
is 



Ja/b, 

A/B, 
A/B 



(Qf 



(Q x ~ v 



= 
= 

= ±d 



In this section, we study the effects of thermal fluctuations 
on the phase diagram for T > 0. We do this using stan- 
dard Weiss mean field theory, taking into account the sym- 



(73) 



metry structure of the phases identified in Sec. Ill The mean 
field method produces self-consistent equations for the order 
parameters, which can be solved, choosing the solution with 
minimal free energy, to obtain the temperature dependence of 
physical quantities. These equations are generally sufficiently 
complicated that, even taking into account symmetry condi- 
tions, only a numerical solution is possible. 

Before presenting the numerical solution, we discuss a few 
aspects of the phase diagram which can be understood analyt- 
ically. Specifically, we consider the instabilities of the param- 
agnetic state on decreasing the temperature from large to small 
values. Several instabilities are possible, which are signaled 
in the mean field calculations by the appearance of a solution 
of a given symmetry S at a temperature T$. Mathematically, 
the temperature T$ is defined by the vanishing of the coeffi- 
cient of the quadratic term in the order parameter associated 
with the symmetry S in the Landau free energy, or equiva- 
lently, the divergence of the mean field susceptibility of the 
S order parameter. The order parameters in question are the 
magnetization vector and quadrupolar tensor at p = and at 
p = Q = 27r(001) (note the four sublattice states should be 
described by several p = 27r(001) wavevectors). In this way 
we obtain four temperatures T m (0), T m (Q), Tq(0), Tq(Q). 
For a given value of the exchange parameters (J, J', V), the 
largest of these temperatures will determine the actual insta- 
bility of the paramagnetic state, and thus the first type of order 
which is encountered upon lowering the temperature. In prin- 
ciple, a first order transition could also occur at a temperature 
higher than this, but the full mean field solution shows that 
this does not occur except in a region where J'/ J, V/J <C 1. 



with q' = 2V3. 



The instability temperatures may determined analytically. 
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One finds 

T m (0) 



6 J' - 2 J V2 J 2 - 20 J J' + 52 J' 2 



7 



Tq(0) = -(J-6J' + 7F), 



15 
7 

180 



J - 2 J' + rcos 



(-3J+ 18 J' +43F). 



T m (Q) 
2g(g) 

Here we defined 

r = V6J 2 -20JJ' + 32J' 2 , 



Q' 



5 - 72y + 204y 2 - 160y 3 



+ i\J{r/Jf - (5 - 72y + 20V - 160y 3 ) ; 
and y = J'/ J. 



(74) 
(75) 
(76) 
(77) 

(78) 
(79) 
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FIG. 2. (Color online) T > phase diagram. In Region I, 
(white) the system transitions directly from the high temperature 
paramagnetic phase to a magnetic state with ordering wavevector 
Q = 2tt(001) at T = T m (Q). In region II (blue), the transition 
from the high temperature paramagnet is instead to a quadrupolar 
phase with ordering wavevector Q, at T — Tq(Q). In region III 
(red), the transition from the paramagnetic is to a q = ferromag- 
netic state at T — T m (0). The dashed curves are the boundaries of 
the ground state phases taken from Fig. [TJ "a, b, c, d" label the low 
temperature phases of each region. J = 1 in the phase diagram. 

In the parameter space depicted in Fig. [T] Tq (p = 0) is al- 
ways smaller than the other three temperatures. Thus there is 
never an instability of the paramagnet to a uniform quadrupo- 
lar ordered state. Comparing the remaining three tempera- 
tures, we find three distinct regions shown in Fig. [2] In region 
I, the highest transition temperature is T m (Q), and magnetic 
order with an enlarged unit cell sets in directly from the para- 
magnetic state. In region II, the highest transition temperature 
is Tq(Q) and two-sublattice quadrupolar (spin nematic) oc- 
curs neighboring the paramagnetic state. In region III, the 
highest transition temperature is T m (Q), and the paramag- 
netic phase undergoes a transition directly to a ferromagnetic 
one. 



On further lowering of temperature, additional phases may 
occur. Full mean-field calculations show that there are in fact 
are ten different patterns of thermal evolution, indicated by the 
the different sub-regions in Fig. [2] We discuss this in further 
detail below. 



A. Region I 



<l/r> 



2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



4=t 



0.1 



0.2 



0.3 



0.4 



0.5 



(a) J' = 0.32 J, V = 0.05 J in I„ 
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(b) J' = 0.25 J, V = in I 6 



<t/C> 




T 

(c) J> = V = 0.1 Jin I c . 

FIG. 3. (Color online) Order parameters plotted for three subregions 

of region I: (a) Sqaure (red) {{Q^'^ - Qf~ y2 })/2, ball (blue) 

{Ql z2 + Qb z2 )/2, diamond (yellow) |(j A + Jb)]/2, and triangle 

(green) |(j A -j B )|/2. (b) Sqaure (red) E 1=a ,b,c,d<Q? z2 >' bal1 ( blue ) 

E l= A.B,c,D<Of ~ y2 )> diamond (yellow) \{f A - j$ + f c - jg)|/4, 
and triangle (blue) \{f h - j% - f c + jg)|/4. (c) Sqaure (red) 

{Q X A~ y2 + Qf~ y2 }/2, ball (blue) (Ql^ + Qf 2 )/2, and diamond 
(yellow) | Oa-Jb) |/2. 
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In region I, the system has a direct transition from the high 
temperature paramagnetic phase to low temperature magneti- 
cally ordered phases specified by the dashed curves in Fig. [2] 
At mean field level, the transitions to FM1 10 phase in region 
I a and AFM100 phase in region I c are found to be first order, 
while a continuous transition to the four-sublattice "A" phase 
is observed in region If, (see Fig. [3]). 



B. Region II 

In region II, there is a broad p = 27r(001) quadrupolar 
phase in the intermediate temperature. Unlike the spin ne- 
matic state in the ground state phase diagram, this intermedi- 
ate temperature quadrupolar phase is actually a biaxial spin 
nematic state in which the quadrupole moment (Q^) has 
three distinct eigenvalues. At mean field level, the transition 
from the PM phase to quadrupolar phase is continuous (see 
Fig.|4]and Fig. [5]). Beyond mean field theory, this transition is 
believed to be in a three dimension 0(3) universality class. 9 

As mentioned previously, the low temperature phases 
(FM110 in region II a , phase "A" in region II C ) can be re- 
garded as further breaking the time reversal symmetry com- 
ing from the quadrupolar phase at intermediate temperature. 
This transition is found to be continuous at mean field level. 
The symmetry breaking associated with this transition can be 
described by several Ising order parameters (uniform and stag- 
gered magnetization). This transition may be continuous be- 
yond mean fieldP 

The transition from quadrupolar phase to AFM100 in re- 
gion lid is found to be strongly first order in mean field theory 
(see Pig. [5j. This is easy to understand from a simple symme- 
try analysis. The AFM100 phase has a uniform orbital config- 
uration and its symmetry is not a subgroup of the quadrupolar 
phase. Hence, it is almost impossible for the transition to be 
continuous. 



C. Region III 

In region III, the intermediate temperature phase is FM1 11. 
The transition from PM to FM111 is found to be continuous 
at mean field leavel. The transition from FM111 to low tem- 
perature phases (FM110 and phase "*") is also found to be 
continuous in our mean field analysis. The transition from 
FM1 1 1 to phase FM1 10 can be simply described by an Ising 
valuable m z . Similarly, the transition from FM111 to phase 
"*" can also be described by an Ising valuable which is the 
staggered magnetization. Thus, these two transitions could be 
continuous beyond mean field. 



D. Magnetic susceptibility 

In this subsection, we discuss the magnetic response at 
T > 0. We focus particularly on the intriguing non-magnetic 
quadrupolar phase of region II in Fig. [2] The magnetic re- 
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T 
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0.5 1.0 1.5 

T 

(b) J' = 0.2 J, V = 0.9 J in II b . 
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T 

(c) J' = 0.05J, V = 1.0J in II C . 
FIG. 4. (Color online) Order parameters plotted in four subregions 

2 _ 2 2 _ 2 

of region II: (a) square (red) (Q x k v — Qg v }/2, ball (blue) 
|(Qa z2 + Qb* 2 }, triangle (green) (j A + j B )|/2, and diamond (yel- 
low) |(j A - j B )|/2. (b) square (red) (Qf - Qt'^)/ 2 and 
ball (blue) \{Q\ z2 + Ql z2 ) (c) square (red) (Qf ~ v * + Qjf^ 2 - 
Q f-v 2 _ Q£ 2 -« 2 }/4, ball (blue) E 1=a ,b,c,d<Q? z2 >/ 4 ' and di " 

2 _ 2 2 _ 2 

amond (yellow) (Q A v + Qg v } /2, upper triangle (green) 
ICJa — Jb)|/2, and down triangle (blue) |(j A +j B )|/2. 



sponse is found to be an important indicator for the quadrupo- 
lar ordering transition. 

The general features, observed in Fig. [7} are as follows. At 
high temperatures, the magnetic susceptibility \ obeys the 
Curie-Weiss law, x -1 ~ A(T - Q cw ), when T > Q cw . 
The Curie- Weiss temperature is readily obtained from a high- 
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FIG. 5. (Color online) Order parameters in region lid with J' = 

0.1J,V = 0.5 J: square (red) (C}f ~ y2 ~ Qf~ y2 }/2, ball (blue) 

\{Q 3 / + Qb 2 }, diamond (yellow) (if - jg + j y - jg>/4. and 
triangle (green) (j" 



if 



temperature series expansion: 



-17J + 66J' 
25 ' 



(80) 



From the mean-field solution, we obtain the susceptibility at 
lower temperature. At the quadrupolar ordering transition, 
a cusp in \ is observed (see Fig. [7). This cusp separates 
the true Curie- Weiss regime of the PM phase from a sec- 
ond Curie-Weiss regime at intermediate temperatures. The 
existence of two Curie- Weiss regimes can be understood as 
due to the remaining magnetic degeneracy of the quadrupolar 
phase. Specifically, the intermediate quadrupolar phase par- 
tially lifts the five-fold spin-orbital degeneracy, giving rise to 
a local doublet which is a time-reversal pair. One should note 
that this pair is not, however, a Kramer's pair. This doublet 
is responsible for the Curie- Weiss behavior at the interme- 
diate temperature regime. Looking in more detail, since the 
spin nematic phase at intermediate temperature has a tetrag- 
onal symmetry, we obtain two different susceptibilities: par- 
allel to the wavevector p = 27r(001) (Xzz) an d normal to it 
(Xxx — Xyy)- As the temperature is lowered further and mag- 
netic order develops, this is reflected in additional features in 
the susceptibility. 



V. DISCUSSION 

In this paper, we introduced and analyzed a spin-orbital 
model to describe the localized electrons in a 4d 2 or 5d 2 con- 
figuration on an fee lattice, in which strong spin-orbit coupling 
and the three-fold degeneracy of the two-electron orbital states 
combine to induce a local effective j = 2 moment. Nearest 
neighbor antiferromagnetic and ferromagnetic exchange inter- 
actions and electric quadrupolar interaction were included in 
the model Hamiltonian. We obtained the ground state and 
finite temperature phase diagrams by Weiss mean field the- 
ory. Seven different ground states (or low temperature) phases 
were found. In a large portion of the parameter space, the sys- 
tem develops a two-sublattice structure of orbital configura- 
tion which is driven by the diagonal orbital-orbital interaction. 
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(c) J = 0.8 J, V = 0.5 J in III C 

FIG. 6. (Color online) Order parameters plotted in three subre- 
gions of region III: (a) uniform magnetization |(j»)| (b) square (red) 

{Qf 2 +Qb z2 }/2, ball (blue) (Qf ~ y2 ~Qf~ y2 )/2, diamond (yel- 
low) |(j A - Jb)/2|, upper triangle (green) |<Ja + j'b )/2|, and down 
triangle (blue) |Ul+jI}/2|. (c) square (red) (Qf -Q'f~ y2 }/2, 

ball (blue) (Q^ 2 + Q% z2 )/2, diamond (yellow) |(j A -j B )/2|, upper 
triangle (green) | (jf + jg ) /2 1 , and down triangle (blue) j (j % + j| ) /2 1 



Most interestingly, a non-magnetic spin nematic ground state 
occurs in the ground phase diagram and extends to a large por- 
tion of the phase diagram at finite temperature. Moreover, we 
find ten different ways for the system to evolve from the high 
temperature paramagnetic phase to the seven low temperature 
phases. 

Our theory has provided numerous preditions for experi- 
ment. For the magnetic ordered phases which break time- 
reversal symmetry, neutron scattering, NMR and/or magne- 
tization measurements can probe the magnetic structure. For 
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FIG. 7. (Color online) Inverse magnetic susceptibility of different 
subregions of region II in Fig. [2j Blue (lower) curve: 1/xaex, red 
(upper) curve: 1/xzz, green (middle) curve: 1/Xvowder- 



the spin nematic phase, similar to the one electron case dis- 
cussed in Ref. [2 the magnetic quadrupole order is expected 
to induce a distortion of the lattice and lower the crystal sym- 
metry. Specifically, the quadrupolar phase corresponds to the 
tetragonal space group PA^/mnm (number 136), , distinct 
from the the cubic space group Fm3m of the high temper- 
ature phase. If this distortion leads to a measurable effect, 
high resolution x-ray scattering should be able to identify the 
spin nematic order. The low temperature spin nematic state 
may also be identified directly from measurements of the or- 
bital state by resonant x-ray scattering or x-ray reflectometry, 
which could be compared with the theoretical wavefunctions 

in Sec. nun 

Now we discuss the specific materials which have been 
studied experimentally to date. We start from P^CaOsOg!^ 
It retains the cubic Fm3m structure down to 17K. The Curie- 
Weiss temperature of this material is — 157K. The magnetic 
moment 1.61/xb is much smaller than the spin only contribu- 
tion for spin S = 1, but close to our prediction 1.25/^b 



based on strong SOC in Sec. IIA The deviation can be un- 
derstood from the effect of hybridization of the Os d orbitals 
with the oxygen p orbitals, which increases the local mag- 
netic moment. 9 Both magnetic susceptibility and specific heat 
measurement find a single antiferromagnetic phase transition 
at Tn = 51K. According to our theory, a single magnetic 
transition correponds to region I in Fig. [2] Since the low tem- 
perature phase of Ba2CaOsC>6 is antiferromagnetic, then the 
AFM100 and "A" phases are consistent with magnetization 
measurements. The enlarged unit cell and detailed orientation 
of magnetic moments predicted here for these phases should 
provide targets for future neutron scattering measurements. 

The structure of La2LiReC»6 was observed to be mono- 
clinic (space group P2i/n)!^This material has a Curie-Weiss 
temperature — 204K, indicating a large antiferromagnetic ex- 
change. The magnetic moment is 1.97/xb, the smallnest of 
which suggests the important of strong SOC. Magnetic sus- 
ceptibility, neutron diffraction and /iSR did not find magnetic 
long range order down to 2K, pointing to a possible quantum 
spin liquid phase in this material. Since this crystal structure 
of this material deviates strongly from a cubic one, our predic- 
tions based on the cubic structure are not applicable. Instead, 
some splitting of the j = 2 manifold should be taken into 
account. Because this is a non-Kramers ion, one could imag- 
ine a trivial non-magnetic ground state at the single ion level. 
However, we note that in Ref fl8l substantial differences were 
observed between zero field and field cooled samples below 
50K, which was argued to be evidence against a single-ion 
singlet. If the local ground state is instead a doublet, an exotic 
ground state could be favored, as proposed for the isostruc- 
tural material L^LiMoOe- 28 More theoretical and experimen- 
tal study of the crystal field splitting and multiplet structure is 
required to further understand the ground state. 

Ba 2 YReC»6 is another material with a cubic crystal struc- 
ture (space group Fm3m).' 8 The Curie-Weiss temperature 
is — 616K, suggestinga a predominant antiferromagnetic ex- 
change. The magnetic moment 1.93/Zb is consistent with 
the picture of strong SOC. The magnetic susceptibility data 
clearly suggests two transitions and shows two Curie regimes. 
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The first transition is at ~ 150K and the second (spin freezing) 
transition is at ~ 50K. The second Curie regime appears at the 
intermediate temperatures between 50K and 150K. Neutron 
diffraction shows the absence of detectable magnetic Bragg 
peaks. /iSR relaxation data observes spin freezing at low tem- 
perature. We may speculate that the spin freezing results from 
disruption by defects of an ordered phase that would otherwise 
occur in an ideal sample. From the existence of two Curie 
regimes, we postulate that Ba2YReOg corresponds to region II 
in the phase diagram Fig. [2] This would identify the interme- 
diate temperature phase as a spin nematic. Interestingly, /iSR 
measured below 100K found evidence for two spin compo- 
nents, which may be consistent with the two-sublattice nature 
of the quadrupolar/spin-nematic phase. Let us consider the 
low temperature phase. Given the very large negative Curie- 
Weiss temperature, the ferromagnetic exchange is likely weak 
in Ba2YReC<6. Hence, comparing with Fig. [T] the natural 
low temperature phases are AFM100, the four-sublattice "A" 
phase or spin nematic phase. Characterization of the type of 
disorder in this material would be helpful in further elucidat- 
ing the physics. It would also be interesting to more directly 
attempt to detect the proposed quadrupolar order experimen- 
tally in the intermediate temperature phase. 

This paper (and the related study in Ref 9) provide a the- 



oretical framework to understand the magnetism and orbital 
physics of this class of materials. Looking to the future, there 
is considerable room for refinement of the theory. It would be 
useful and interesting to include the Jahn-Teller effect, and to 
develop some microscopic understanding (perhaps from ab- 
initio calculations) of the crystal field splittings in non-cubic 
materials. Within the present model, more studies of ther- 
mal and quantum fluctuations beyond mean-field and spin- 
wave approaches would be desirable. An understanding of 
the types of disorder in these materials and their effects on the 
magnetism is also needed. Hopefully continued pursuit and 
refinement of measurements on this interesting class of mate- 
rials will motivate further theoretical work on these and other 
points. 
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